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Abstract 


This  report  summarizes  research  conducted  during  the  past  three  and  a  half 
years  aimed  at  developing  and  testing  a  turbulence/transition  model  applicable 
to  high-speed  turbulent  hows.  The  first  two  years  of  the  project  focused  on 
fully  turbulent  flows,  while  emphasis  shifted  to  boundary-layer  development  in 
the  transition  region  during  the  final  year  and  a  half.  This  report  includes  a 
brief  summary  of  research  accomplished  during  the  first,  three  years  and  cites 
publications  that,  describe  research  results  in  greater  detail.  The  main  body 
of  this  report  summarizes  research  conducted  during  the  final  six  months  of 
the  period  of  performance.  The  primary  results  of  the  last  six  months  of  the 
project  are  elimination  of  the  k-u  model’s  sensitivity  to  the  freestream  value  of 
u  and  development  of  a  method  for  triggering  transition  at  a  specified  location, 
independent  of  the  freestream  turbulence  level. 
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Chapter  1 


Introduction 


The  purpose  of  this  research  project  for  the  past  three  and  a  half  years  has  been 
to  develop  analytical  and  computational  tools  suitable  for  predicting  properties 
of  hypersonic  flows,  including  both  turbulent  and  transitional  flow  regimes.  Our 
research  efforts  have  generated  a  total  of  eight  publications,  including  this  final 
report.  The  previous  publications  are  as  follows. 

1.  Wilcox,  D.  C.,  “Hypersonic  Turbulence  Modeling  Without  the  Epsilon 
Equation,”  1th  National  Aero-Space  Plane  Technology  Symposium,  Ses¬ 
sion  VII,  Paper  No.  27  (October  1989). 

2.  Wilcox,  D.  C.,  “A  Half  Century  Historical  Review  of  the  k-uj  Model,” 
AIAA  Paper  91-0615  (January  1991). 

3.  Wilcox,  1).  C.,  “Theoretical  Study  of  'IXirbulent  Mixing  Between  Hyper¬ 
sonic  Streams,”  DCW  Industries  Report  No.  DCW-R-35-01  (May  1991). 

4.  Wilcox,  I).  C.,  “Progress  in  Hypersonic  'IWbuIence  Modeling,”  AIAA  Pa¬ 
per  91-1785  (June  1991). 

5.  Wilcox,  1).  C.,  “Progress  in  Developing  a  Transition  Model  for  High-Speed 
Flows,”  DCW  Industries  Report  No.  DCW-R-35-02  (October  1992). 

6.  Wilcox,  D.  C.,  “The  Remarkable  Ability  of  Turbulence  Model  Equations  to 
Describe  Transition,”  f  ifth  Symposium  on  Numerical  and  Physical  Aspects 
of  Actvdynamic  Flows,  California  State  University,  Long  Beach,  California 
(13-15  January  1992)  [accepted  for  publication  in  the  AIAA  Journal]. 

7.  Wilcox,  D.  C.,  “Dilatation-Dissipation  Corrections  for  Advanced  Turbu¬ 
lence  Models,”  AIAA  Journal,  Vol.  30,  No.  11,  pp.  2639-2646  (November 
1992). 


Publications  1,  3,  4  and  7  deal  primarily  with  fully-turbulent,  high-speed 
flows,  and  represent  research  results  for  the  first  two  years  of  the  project.  One 
of  the  primary  accomplishments  was  discovery  that  separating  dissipation  into 
solenoidal  and  dilatation  contributions  improves  predictive  accuracy  for  the 
compressible  mixing  layer  at  the  expense  of  a  loss  in  predictive  accuracy  for 
boundary  layers.  We  identified  the  cause  of  this  problem  and  postulated  a 
dilatation-dissipation  correction  that  is  accurate  for  the  mixing  layer  and  the 
boundary  layer.  A  second  accomplishment  was  identification  that  the  k-u)  model 
faithfully  reproduces  the  compressible  law  of  the  wall  for  high-speed  flows  while 
the  Jfc-<  model  does  not.  Perhaps  the  most  important  accomplishment  was  the 
demonstration,  in  Publications  4  and  7  that  dilatation  dissipation  does  little  to 
improve  two-equation-model  predictions  for  shock-separated  flows.  Rather,  the 
research  verifies  that  using  a  second-order  closure  model  greatly  improves  pre¬ 
dicted  size  of  the  separation  region,  although  other  flow  properties  are  not  accu¬ 
rately  predicted  (e.g.,  velocity  profiles  downstream  of  the  reattachment  point). 

Publications  2,  4  and  6  introduce  iow-Reynolds-number  corrections  for  the 
k-u)  model  that  permit  the  model  to  accurately  describe  near-wall  properties 
of  high-speed  houndary  layers.  The  modifications  presented  in  the  first  two 
publications  ultimately  have  been  superseded  by  those  postulated  in  Publication 
6.  The  latter  Publication  represents  much  of  the  basic  research  conducted  during 
the  third  year  of  the  project. 

Publications  5  and  6  focus  on  application  of  the  low-Reynolds-uumber  k-u) 
model  to  a  collection  of  transitional  flows.  The  publications  show  that  the 
model  accurately  predicts  properties  in  the  transitional  region  provided  tran¬ 
sition  Reynolds  number  is  relatively  small.  By  contrast,  predicted  transition 
width  is  generally  smaller  than  measured  for  transition  at  high  Reynolds  num¬ 
bers.  Additionally,  computed  results  show  a  stronger-than-desired  sensitivity 
to  the  freestream  value  of  w,  and  the  need  for  a  simple  method  of  triggering 
transition  at  a  specified  location. 

During  the  final  six  months  of  the  project,  we  have  eliminated  two  of  the 
three  outstanding  problems  with  the  model  as  it  applies  to  transitional  bound¬ 
ary  layers.  Specifically,  we  have  devised  a  modification  that  eliminates  solution 
sensitivity  to  the  freestream  value  of  w,  and  we  have  devised  a  straightfor¬ 
ward  method  for  triggering  transition  at  a  specified  location.  However,  limited 
funding  has  precluded  solving  the  problem  of  the  too-abrupt  transition  at  high 
Reynolds  numbers. 

As  summarized  in  Publication  5,  the  overall  objective  of  our  approach  to 
the  transition  problem  is  to  use  the  Wilcox1  k-u>  turbulence  model  as  the  foun¬ 
dation  for  studying  and  modeling  the  transitional  flow  region.  Consistent  with 
the  needs  of  NASA,  the  transition  point  is  assumed  to  be  known  a  priori.  Com¬ 
pulations  can  thus  be  initiated  at  the  known  transition  location  and  continued 
downstream  through  the  transitional  flow  region  and  well  into  the  fully  turbu¬ 
lent  region.  To  develop  the  model,  we  have  followed  a  sequence  of  interrelated 
steps. 
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1.  In  Publications  2,  4  and  6,  we  drew  from  the  extensive  work  done  by 
Wilcox  et  al.2-7  to  help  formulate  a  low-Reynolds-number  version  of  the 
k-ui  model. 

2.  In  Publications  4  and  6,  we  simulated  several  of  the  flows  that  have  been 
done  with  Direct  Numerical  Simulation  (DNS)  methods.  The  results  of 
Mansour,  Kim  and  Moin8  were  computed  for  two  different  Reynolds  num¬ 
bers.  Comparing  model  predictions  with  the  DNS  results  helped  greatly 
in  developing  the  model. 

3.  In  Publication  5,  we  tested  the  model  against  10  two-dimensional,  fully 
turbulent  boundary  layers,  including  both  incompressible  and  compress¬ 
ible  cases.  We  have  repeated  5  of  these  cases  in  this  report. 

4.  In  Publication  5,  we  tested  the  model  against  all  of  the  two-dimensional 
cases  presented  by  Singer,  Dinavahi,  and  Iyer.9  We  have  repeated  all  of 
these  cases  in  this  report. 

5.  We  have  analyzed  the  compressible  log  layer  to  explain  why  adding  a  cross- 
diffusion  term  to  the  u>  equation  can  have  an  adverse  effect  on  predicted 
skin  friction. 

As  noted  above,  addressing  Items  1  and  2  has  produced  a  low-Reynolds- 
uumber  version  of  the  k-u>  model.  Publication  6  describes  the  model  and 
presents  applications  to  fully  turbulent  channel  and  pipe  flow,  and  for  a  tran¬ 
sitional,  incompressible  flat-plate  boundary  layer.  In  this  report  we  propose  a 
modified  version  of  the  model  that  eliminates  freestream  boundary-condition 
sensitivity.  This  report  applies  the  model  to  5  turbulent  boundary-layer  test 
cases,  two  free  shear  flow  cases  and  to  more  than  20  transitional  flows  for  which 
experimental  data  are  available. 

Chapter  2  summarizes  the  revised  low-Reynolds-number  version  of  the  k-ui 
model.  Chapters  3,  4  and  5  present  results  of  our  applications.  Chapter  6 
summarizes  results  of  the  research.  The  effect  of  a  model  revision  on  the  com¬ 
pressible  law  of  the  wall  is  discussed  in  Appendix  A. 
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Chapter  2 

Revised  Low  Reynolds 
Number  k-a;  Model 


'Hie  revised  low-Roynolds-number  version  of  the  k-ui  model  is  written  in  terms 
of  Favre10  mass  averaged  variables.  The  equations  of  motion  are  as  follows. 


Op  ,  Or,, 


I  w = £-  [*<’*<  -*+(/*+ <2 
1 (,,t) + £-  (<>«.*) = [<" + «-)^]  <2 


$  .  .  .  uj  T  dui  „  2  P  Ok  Oui 

S(/«)+S-(Wu.)=aru_-^ 

#  r.  , 


In  Equations  (2.1)-(2.5),  <  and  ar*  denote  time  and  position  vector;  p  and 
p  are  density  and  pressure;  «,  is  mass-averaged  velocity  vector;  E  and  H  are 
total  energy  and  total  enthalpy;  and  are  the  total  stress  tensor  and  the 
total  heat  flux  vector;  t?  is  the  Reynolds  stress  tensor;  k  and  u  are  turbulence 


kinetic  energy  and  specific  dissipation  rate;  and,  a,  /J ,  fim ,  a,  a’ ,  04  are  closure 
coefficients.  The  following  constitutive  relations  are  needed  to  close  the  system. 


*  fj  T 

PT  =  Ot  '""l"  ,  Tijt  Tjj  +  T^j 

Ui 

(2.6) 

T y  =  2p  (5” ~  zdrk6ij) 

(2.7) 

rJ  =  2pT  (s’y  -  _ 

(2.8) 

<* 

II 

l 

+ 

51-b 

ns* 

(2.9) 

=  t  +  +  k,  H  =  h  +  -u,u,  +  k 

(2.10) 

1  /  \ 

^  “  2 

(2.11) 

In  the  constitutive  relations,  we  have  introduced  the  molecular  stress  tensor, 
t/j,  mean  strain-rate  tensor,  $y,  molecular  viscosity,  p,  eddy  viscosity,  /it. 
enthalpy,  h,  laminar  Prandtl  number,  Pr*,,  turbulent  Prandtl  number,  Prr, 
internal  energy,  e,  and  an  additional  closure  coefficient,  a* .  Finally,  to  complete 
closure  of  the  system,  we  specify  the  values  of  the  closure  coefficients  as  follows. 


*  a*  +  ReT/Rk  1  Oo  +  Rer/Ru  ,  ,x_i 

°  “  1  +  Rer/Rt  ’  °  ~  2  '  1  +  ReT/R,  ' K  ’ 


(2.12) 


b/\8  +  (ReT/R0)* 
P  100  '  1  +  (RcT/Rp)* 


<*d  - 


0, 


dk  du 

dij  dxj 


<  0 


3_  dk  duj 

10  ’  dxj  dxj  > 


^  -  0/3,  a0  =  1/10 


(2.13) 

(2.14) 

(2.15) 

(2.16) 


Rp  =  6,  Rk  =  8,  Rw  =  11/5 

The  quantity  Rer  is  the  turbulence  Reynolds  number  defined  by 


(2.17) 
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(2.18) 


The  primary  difference  between  this  model  and  the  model  introduced  by 
Wilcox11  is  the  introduction  of  the  term  in  Equation  (2.5)  proportional  to  ad. 
This  term  is  similar  to  Menter’s12  cross-difTusion  term.  The  effect  of  this  term  is 
to  replace  the  entrainment  velocity,  v,  in  the  u>  equation  by  (v  —  <Tdu>~xdk/dy). 
Since  k  decreases  approaching  the  shear  layer  edge  (assuming  rrd  >  0),  the 
net  effect  is  to  make  the  effective  entrainment  velocity  positive  (or  at  least  less 
negative).  As  a  result,  u  diffuses  from  the  turbulent  region  into  the  nonturbulent 
region,  which  is  the  opposite  of  what  happens  with  the  k-ui  model.  Thus,  the 
freestream  value  of  u  has  no  effect  on  the  solution. 

As  shown  by  Wilcox13,  it  is  important  to  suppress  the  cross-diffusion  term 
close  to  solid  boundaries  for  wall-bounded  flows.  The  prescription  in  Equa¬ 
tion  (2.15)  causes  ad  to  vanish  near  a  solid  boundary  since  k  increases  and  u> 
decreases  in  the  viscous  sublayer.  This  modification  to  the  u>  equation  eliminates 
the  model’s  sensitivity  to  the  freestream  value  of  u>.  The  values  of  several  clo¬ 
sure  coefficients  must  be  modified  to  achieve  optimum  results  for  both  boundary 
layers  and  free  shear  layers.  Specifically,  or  and  it*  must  assume  values  larger 
than  used  in  the  Standard  k-ui  model.  This,  in  turn,  requires  changing  or  and 
R u,  to  maintain  a  satisfactory  law  of  the  wall  in  the  turbulent  boundary  layer. 
The  value  quoted  for  a  guarantees  that  the  Karman  constant  is  0.41  while  the 
value  of  Ru  yields  an  additive  constant  in  the  law  of  the  wall  of  5.0.  In  summary, 
Table  2.1  shows  how  the  revised  model’s  closure  coefficients  dhTer  from  those  of 
the  Wilcox11  model. 


Table  2.1:  Closure  Coefficient  Differences 


Coefficient 

New  Model 

Wilcox11  Model 

<T 

3/5 

1/2 

IT* 

1 

1/2 

<fd 

Equation  (2.15) 

0 

a 

1/2 

5/9 

Rw 

11/5 

27/10 

Chapter  3 

Turbulent  Boundary-Layer 
Applications 


A  key  objective  of  research  during  the  final  year  and  a  half  of  this  project  has 
been  to  describe  boundary-layer  development  through  transition  from  laminar 
to  turbulent  flow.  Of  course,  it  is  important  that  we  provide  an  accurate  de¬ 
scription  in  the  turbulent  region  inuncdiateiy  following  transition,  i.e.,  we  insist 
that  our  model  approach  the  proper  limiting  state  of  the  boundary  layer.  Con¬ 
sequently,  since  the  k-w  model  without  viscous  corrections  generally  is  very 
accurate  for  turbulent  boundary  layers,  a  round  of  tests  is  in  order  to  make  sure 
model  predictions  are  not  adversely  affected  by  the  viscous  corrections.  We 
have  performed  five  boundary  layer  computations  including  effects  of  adverse 
and  favorable  pressure  gradient,  and  for  Mach  numbers  as  high  as  10.  Table  3.1 
lists  the  five  cases. 

Table  3.1:  Turbulent  Boundary  Layer  Test  Cases 


Description 

Vp 

Data  Source 

Incompressible  Flat  Plate 

None 

Wieghardt-TillmanM 

Samuel-Joubert 

Adverse 

Samuel-Joubert15 

Mach  2.244  Flat  Plate 

No::" 

Shutts16 

Mach  4.544  Flat  Plate 

None 

Coles16 

Mach  10.31  Flat  Plate 

None 

Watson16 

Incompressible  Flat  Plate.  Figure  3.1  compares  computed  and  measured 
skin  friction,  c/,  and  sublayer-scaled  velocity,  u+.  As  shown,  the  only  significant 
difference  between  predictions  with  and  without  viscous  corrections  is  in  the 
transition  point.  The  unmodified  model  undergoes  transition  much  earlier  than 
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the  model  with  viscous  corrections.  Skin  friction  at  the  final  station  is  2.51  ■  10-3 
with  viscous  modifications  compared  to  2.49  •  10-3  without. 

Samuel- Jouhert .  This  incompressible  adverse  pressure  gradient  boundary 
layer  was  supposed  to  be  a  simple  application  in  Stanford  Olympics  II  17  On 
the  contrary,  it  proved  to  be  very  difficult  for  all  models  and  has  become  a  key 
test  case  for  how  well  a  turbulence  model  predicts  effects  of  adverse  pressure 
gradient,  figure  3.2  shows  that  both  high-  and  low-Reynolds-nuinber  versions 
of  the  model  nearly  duplicate  measured  skin  friction  and  velocity  profile  at  tin- 
last  station.  At  the  final  station,  cj  =  1.35  ■  ID-3  with  viscous  corrections  and 
1.25  •  1(T3  without. 

Compressible  Flat  Plates.  As  shown  in  Figures  ,.3,  3.4  and  3.5,  the 
viscous  corrections  have  virtually  no  effect  on  skin  friction  and  velocity  proliles 
in  the  fully  turbulent  region  for  these  three  applications.  For  all  three  cases, 
skin  friction  differs  by  less  than  one  percent  at  the  final  station. 

In  summary,  for  all  5  cases  considered,  differences  between  computed  How 
properties  with  and  without  viscous  corrections  are  less  than  M%,  and  are  gen¬ 
erally  less  than  2%.  Additionally,  although  not  shown  here,  we  have  computed 
(under  separate  funding)  14  more  incompressible  boundary  layers,  3  with  fa¬ 
vorable  pressure  gradient  and  11  with  adverse  pressure  gradient.  Along  with 
the  flat-plate  boundary  layer  and  the  Samuel-Joubert  case,  these  are  the  16 
boundary  layers  analyzed  by  Wilcox18.  As  with  the  cases  discussed  above,  the 
low-Reynolds-nuinber  fb-ui  model  is  as  accurate  as  the  basic  model  without  vis¬ 
cous  modifications.  Hence,  the  viscous  corrections  leave  the  best  features  of 
the  ir-w  model  intact,  i.e.,  the  model  still  ;»ccurately  predicts  effects  of  pressure 
gradient  and  compressibility  up  to  Mach  5. 

While  the  cross-diffusion  term  in  the  u;  equation  greatly  reduces  the  model's 
sensitivity  to  the  freestream  value  of  w,  it  does  not  entirely  remove  it.  A  larger 
value  of  "i  is  needed  to  completely  remove  the  sensitivity.  However,  as  shown 
in  Appendix  A,  using  a  larger  value  for  04  yields  an  inaccurate  prediction  for 
the  compressible  law  of  the  wall. 
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Figure  3.1:  Incompressible  (lal- plate  boundary  layer;  o  Wieghardt-Tillman 
- k-w  with  viscous  corrections; - k-w  without  viscous  corrections. 


Figure  3.2:  Incompressible  boundary  layer  in  an  increasingly  adverse  pressure 

gradient;  o  Samuel-Joubert; - k-ui  with  viscous  corrections; - k- u>  without 

viscous  corrections. 
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Figure  3.3:  Mach  2. ‘244  adiabatic-wall  flat-plate  boundary  layer;  o  Shults 
- k-w  with  viscous  corrections; - k-uj  without  viscous  corrections. 
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Chapter  4 

Free  Shear  Flow 
Applications 


While  k-u>  model  boundary-layer  computations  display  a  mild  sensitivity  to  the 
freestream  value  of  ui,  the  sensitivity  is  much  stronger  for  free  shear  flows.  As 
shown  by  Wilcox13,  the  spreading  rate  for  the  incompressible  mixing  layer  varies 
between  0.103  and  0.141;  the  measured  spreading  rate  is  0.115.  Similarly,  for 
the  far  wake,  computed  spreading  rate  varies  from  0.301  to  0.500,  as  compared 
to  a  measured  value  of  0.365. 

The  new  model  almost  completely  eliminates  this  problem.  Figure  4.1  com¬ 
pares  computed  and  measured  mixing-layer  velocity  profiles  for  two  values  of 
the  freestream  ui  that  differ  by  several  orders  of  magnitude.  The  computed 
spreading  rate  varies  between  0.105  and  0.107. 

Similarly,  Figure  4.2  compares  computed  and  measured  velocity  profiles  for 
the  far  wake.  For  a  wide  range  of  freestream  values  of  w,  the  spreading  rate 
varies  between  0.295  and  0.318. 

We  have  attempted  computations  for  the  plane  and  round  jets,  but  have 
encountered  numerical  difficulties  in  obtaining  converged  solutions.  This  illus¬ 
trates  a  key  problem  attending  use  of  the  cross- diffusion  term.  In  general,  the 
term  greatly  increases  the  stiffness  of  the  model  equations  by  altering  the  effec¬ 
tive  time  scale  (w- 1 )  near  a  turbulent/noulurbuient  interface.  As  a  consequence, 
the  convergence  rate  of  the  numerical  algorithm  tends  to  be  reduced. 
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Figure  4.1:  Comparison  of  computed  and  measured  incompressible  mixing-layer 
velocity  profiles;  o  Liepmann-Laufer; - <*>«,  — ►  0; - u/oo  finite. 
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Figure  4.2:  Comparison  of  computed  and  measured  incompressible  far-wake 
velocity  profiles;  o  Fage-Falkner; - — *  0; - ui^  finite. 
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Chapter  5 

Transitional 

Boundary-Layer 

Applications 


5.1  Triggering  Transition 

The  question  of  how  to  trigger  transition  for  a  given  application  raises  an  impor¬ 
tant  point  about  using  turbulence-model  equations  for  transitional  Hows.  We 
must  keep  in  mind  that  transition  is  a  complicated  phenomenon.  lYansition 
is  triggered  by  a  disturbance  in  a  boundary  layer  only  if  the  frequency  of  the 
disturbance  falls  in  a  specific  band.  Reynolds  averaging  has  masked  all  spectral 
effects,  and  all  the  model  can  represent  with  k  and  u>  is  the  intensity  of  the 
disturbance  and  an  average  frequency.  Hence,  it  is  possible  for  the  turbulence 
model  to  predict  transition  when  it  shouldn’t  occur.  If  we  choose  to  trigger 
transition  via  the  presence  of  a  freestream  disturbance,  using  turbulence-model 
equations  is  sensible  only  if  the  triggering  disturbance  is  broad  band,  i.e.,  con¬ 
tains  all  frequencies. 

In  devising  the  model,  the  closure  coefficients  a*  and  a0  [Equation  (2.16)] 
have  been  chosen  to  match  the  minimum  critical  Reynolds  number  for  the  Bla- 
sius  boundary  layer.  Consequently,  the  model  can  be  expected  to  yield  sensi¬ 
ble  transition  predictions  for  constant-pressure,  incompressible  boundary  layers. 
Unfortunately,  the  model  equations  must  be  recalibrated  (i.e.,  a*  and  a0  must 
be  adjusted)  to  accommodate  each  new  complicating  effect  [sec  Wilcox  ct  al.2-7]. 

In  this  project,  the  goal  has  been  to  assume  the  transition  point  is  given  and 
to  use  the  model  equations  through  the  transitional  region  and  into  the  turbulent 
part  of  the  flow.  To  date,  the  freestream  turbulence  intensity,  proportional  to 
the  freestream  value  of  k,  has  been  adjusted  to  match  the  measured  transition 
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point.  This  is  satisfactory  when  the  transition  point  occurs  at  a  large  Reynolds 
number,  which  requires  *<»  to  be  small  relative  to  !/£,.  However,  for  high-speed 
flows  in  which  transition  occurs  at  a  relatively  small  Reynolds  number,  we  have 
found  that  unreasonably  large  values  of  koa  are  needed  to  trigger  transition,  so 
large  as  to  affect  the  total  energy  in  the  freestream  in  a  physically  unrealistic 
manner.  Thus,  a  new  method  for  triggering  transition  is  needed. 

To  devise  an  alternative  method,  we  can  take  advantage  of  a  unique  feature 
of  the  k-ui  model.  Specifically,  by  using  a  finite  value  for  u >  at  the  surface,  we 
can  simulate  surface  roughness  with  the  model.  For  fully  turbulent  boundary 
layers,  Wilcox1  shows  that 


Ui  = 


at  y  =  0 


(5.1) 


where  uT  is  friction  velocity  and  Sr  is  a  dimensionless  function  of  the  surface 
roughness  height,  Icr,  defined  by  (with  k^  =  uT*/?/i /): 


Sr=  1 


(50/*+ )2, 

100/*+, 


*+<25 
*«  >  25 


(5.2) 


Since  increasing  the  surface  roughness  height  corresponds  to  decreasing  the  sur¬ 
face  value  of  w  (and  thus  the  dissipation  in  the  *  equation),  the  model  predicts 
that  roughness  will  have  a  destabilizing  effect.  This  is  consistent  with  measure¬ 
ments,  and  patches  of  surface  roughness  are  often  used  to  trigger  transition  in 
experiments.  Thus,  a  possible  way  to  trigger  transition  with  the  model  equa¬ 
tions  is  to  numerically  simulate  a  roughness  strip  via  Equations  (5.1)  and  (5.2). 

We  have  run  more  than  20  two-dimensional  transitional  boundary  layer  cases 
to  test  this  idea;  results  of  the  computations  Eire  given  in  the  next  section.  We 
have  been  able  to  trigger  transition  at  the  desired  location  for  all  of  the  rases 
considered  using  a  roughness  strip  with  *j?  and  the  streamwise  extent  of  the 
strip,  As,  given  by  the  following  correlations. 


^ -  =  0.0lby/R^  (5.4) 

°t 

The  quantities  bt  and  He,t  are  the  boundary-layer  thickness  and  transition 
Reynolds  number  based  on  arc  length. 


5.2  Applications 

In  order  to  test  the  new  transition-triggering  method,  we  have  computed  ail  of 
the  two-dimensional  transitional  boundary  layer  cases  considered  by  Singer,  et 
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al.9  Our  applications  also  include  a  Mach  20  re-entry  case.  Table  5.1  summarizes 
the  cases  we  have  done. 

Table  5.1:  Transitional  Boundary  Layer  Test  Cases 


Flow 

Description 

Vp 

Data  Source 

Incompressible  Flat  Plate 

None 

Schubauer-Klebanoff19 

Favorable  Vp 

Favorable 

Blair-Werle20 

Supersonic  Cone  Flow 

None 

Fisher-  Dougherty21 

Supersonic  Cone  Flow 

None 

Chen,  et  al.22 

Freestream  Turbulence 

None 

Schubauer-Skramstad23 

Freestream  Ihrbulence 

None 

Blair24 

Prolate  Spheroid 

Favorable 

Meier- Kreplin-Ming25 

Concave  Surface 

None 

S  weari  ngen-  Black  welder26 

6 

Mach  20,  5°  Cone 

None 

Howard27 

In  all  cases,  computation  begins  at  the  plate  leading  edge,  and  the  turbu¬ 
lence  kinetic  energy  is  initially  set  to  an  extremely  small  value,  viz.,  10-15i/£,, 
throughout  the  boundary  layer.  This  value  is  too  small  to  trigger  transition 
naturally.  The  initial  ui  profile  is  given  by  the  exact  laminar-flow  solution  to  the 
model  equations.11 

In  addition  to  testing  the  transition  triggering  method  described  above,  we 
have  tested  the  effect  of  the  cross-diffusion  term  in  Equation  (2.5).  A  disap¬ 
pointing  feature  of  our  previous  transition  predictions  was  a  mild  sensitivity  of 
transition  location  to  the  initial  u  profile.  The  sensitivity  was  caused  by  the 
small,  but  noticeable,  effect  of  the  freestream  value  of  w  on  the  laminar-flow  so¬ 
lution  for  w.  The  primary  reason  for  introducing  the  cross-diffusion  term  was  to 
eliminate  this  sensitivity  for  transition  predictions.  As  shown  in  the  preceding 
chapters,  the  cross-diffusion  term  eliminates  most  of  the  sensitivity  for  turbu¬ 
lent  shear  flows.  Numerical  experimentation  for  several  of  the  transitional  cases 
listed  in  Table  5.1  verifies  that  this  is  true  for  transitional  flows  also.  That  is, 
our  computations  verify  that  transition  predictions  are  completely  insensitive  to 
the  initial  u  profile,  even  when  transition  is  triggered  by  freestream  turbulence. 

Incompressible  Flat  Plate.  The  first  case  is  Flow  1  from  Singer,  et  al.9 
This  is  an  incompressible  flat-plate  boundary  layer  that  undergoes  transition 
at  a  plate  length  between  1.6  m.  and  1.8  m.  According  to  Equations  (5.3)  and 
(5.4),  this  flow  requires  a  roughness  strip  with  kn/6t  —  3.0  and  A s/6t  =  25. 
Numerical  experimentation  shows  that  a  shorter  transition  strip,  i.e.,  A  s/6,  =  7 
is  sufficient  to  trigger  transition  at  the  desired  location  for  this  case. 

Figure  5.1  compares  computed  and  measured  skin  friction  throughout  the 
transition  region  for  the  model  with  and  without  the  cross-diffusion  term.  In 
both  cases,  computed  and  measured  cj  differ  by  less  than  16%  of  the  peak  skin 
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Figure  5.1:  Transitional  incompressible  flat-plate  boundary  layer; - New  k-ui 

model;  -  -  -  k-ui  model  without  cross  diffusion. 


friction.  Note  that  with  cross-diffusion  included,  the  distance  between  minimum 
and  maximum  skin  friction  is  a  bit  less  than  that  without  cross  diffusion. 

Favorable  Pressure  Gradient  Boundary  Layers.  The  next  applications 
are  for  incompressible  boundary  layers  in  a  favorable  pressure  gradient.  The 
boundary  layers  considered  correspond  to  Flow  2/Cases  1  and  2  of  Singer,  et  al. 
In  addition  to  having  adverse  pressure  gradient,  the  surface  is  cooled.  Figure  5.2 
compares  computed  and  measured  Stanton  number  for  the  two  rases.  Case  2 
has  a  stronger  favorable  pressure  gradient  than  Case  1.  The  dimensions  of  the 
roughness  strip  required  to  match  the  measured  transition  point  are  (kn/6t, 
As/fit)  =  (10,  10)  and  (4,  9),  respectively,  for  Cases  1  and  2.  Equations  (5.3) 
and  (5.4)  indicate  (knf&i,  As/6t )  =  (8.5,  8.7)  and  (6.0,  12.5),  respectively. 
Again,  we  see  that  the  transition  width  is  reduced  a  bit  when  the  cross-diffusion 
term  is  included. 

Supersonic  Cone  Flow.  For  this  application  we  consider  Flow  3  of  the 
Singer,  et  al.  study.  Cases  1-4  focus  on  flow  past  a  5°  half  angle  cone  with 
Mach  numbers  ranging  from  1.16  to  1.86,  corresponding  to  measurements  of 
Fisher  and  Dougherty.21  Cases  5-7  are  for  Mach  3.36  flow  past  a  5°  half-angle 
cone,  with  an  adiabatic  surface  (see  Chen,  et  al.22).  Table  5.2  summarizes  the 
transition  Reynolds  number,  KeXi,  and  Reynolds  number  based  on  transition 
width,  for  all  7  cases,  for  the  new  and  old  low-Reynolds-number  k-ui 

models. 
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Case  2:  Strongly  favorable  pressure  gradient 


Figure  5.2:  Transitional  incompressible  boundary-layer  flow  with  favorable  pres¬ 
sure  gradient  and  surface  cooling; - New  k- w  model; - k-u>  model  without 

cross  diffusion. 
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Table  5.2:  Results  for  Fisher-Dougherty  and  Chen  et  al.  Test  Cases 


Case 

Mach  Number  He 

£( 

(  Rf  Aft  )nrui 

{He  Ax,)old 

(  )f!TV 

1 

1.16 

7.13 

10* 

0.45 

10e 

0.89 

W~ 

0.69 

To5- 

2 

5.70 

106 

0.50 

106 

0.84 

106 

0.97- 

10® 

3 

1.55 

7.90 

106 

0.72 

106 

1.08 

106 

1.24 

10® 

4 

1.86 

7.49 

106 

0.78 

106 

1.07 

106 

1.56- 

106 

5 

3.36 

6.92 

106 

1.01 

106 

1.63 

106 

— 

6 

3.36 

9.51 

106 

1.19 

106 

1.79 

106 

— 

7 

3.36 

1.07 

107 

1.30- 

10s 

1.82- 

106 

— 

Freestream  Turbulence  Effects.  We  turn  now  to  effects  of  frcestream 
turbulence,  i.e.,  Flow  \  of  the  Singer,  et  al.  study.  All  of  the  cases  considered  are 
for  incompressible  boundary  layers.  Table  5.3  shows  that  computed  transition 
w>dth  for  Cases  1-5  is  consistently  smaller  than  measured,  and  noticeably  smaller 
than  the  model  without  cross  diffusion. 


Table  5.3:  Results  for  Schubauer-Skramstad  Test  Cases 


PP31 

mm 

Rext 

( Re &Tl  )oid 

(^Ar,  )f/p 

mm 

.042% 

2.85  10e 

0.56  10* 

0.88  •  10* 

1 .00  10® 

mm 

.100% 

2.75  •  10® 

0.52-  10® 

0.86  - 10® 

1.20  10® 

Kfl 

.200% 

2.20  •  106 

0.48  •  10® 

0.78  10® 

1.50-10® 

m 

.260% 

1.80  106 

0.45  •  10® 

0.69-  10® 

1.40  10® 

mm 

.340%, 

1.40  10® 

0.36  10® 

0.58-  10® 

1.25-10® 

For  Cases  6,  7  and  8,  Figure  5.3  shows  that  predicted  peak  skin  friction  is 
about  10%-15%  lower  than  measured.  The  results  for  these  three  cases  are  very 
similar  to  those  obtained  without  cross  diffusion. 

Incompressible  Flow  Past  a  Spheroid.  This  case  is  Flow  5  from  Singer, 
et  al.  The  flow  examines  transition  of  a  boundary  layer  over  a  prolate  spheroid 
at  zero  angle  of  attack,  with  transition  triggered  by  a  roughness  strip.  Fx- 
perimental  data  have  been  provided  by  Meier,  Kreplin  and  Miug.25  Table  5.4 
summarizes  the  four  cases. 

Figures  5.4  through  5.7  compare  computed  and  measured  skin  friction  (based 
on  local  boundary-layer  edge  velocity)  as  a  function  of  arc  length  along  the 
spheroid.  As  shown,  computed  and  measured  skin  friction  are  closest  when 
transition  is  induced  by  the  roughness  strip.  Without  the  roughness  strip,  the 
predicted  transition  occurs  more  abruptly  than  measured. 
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Table  5.4:  Meier-Kreplin-Ming  Test  Cases 


Case 

Description 

Er/6, 

As/6, 

1 

Uoo  =  20  m/scc,  no  roughness  strip 

2.0 

6.0 

2 

Uoo  --  20  m/sec,  with  roughness  strip 

4.0 

8.0 

3 

Uoo  —  30  m/sec,  no  roughness  strip 

2.0 

1.6 

4 

f/oo  =  30  in/sec,  with  roughness  strip 

2.0 

15.0 

Coucave  Surface  Boundary  Layer.  In  this  application  we  consider  in¬ 
compressible  flow  over  a  concave  surface.  We  present  results  in  Figure  5.8.  As  in 
all  previous  applications,  including  cross  diffusion  causes  a  noticeable  reduction 
in  transition  width.  According  to  Equations  (5.3)  and  (5.4),  the  transition  strip 
should  have  ifc/,/6,  =  3.5  and  As/6,  =  21.  Numerical  experimentation  shows 
that  somewhat  smaller  values,  E„/6,  =  3.0  and  As/6,  =  12,  can  be  used. 

Mach  20,  5°  Cone.  Figure  5.9  compares  computed  and  measured  Stan¬ 
ton  number,  St,  for  Mach  20  flow27  past  a  5°  half-angle  cone.  As  shown,  the 
computed  Stanton  number  increases  more  abruptly  than  measured,  although 
less  abruptly  with  cross  diffusion  included.  This  is  the  only  case  for  which 
transition  width  increases  when  cross  diffusion  is  included.  Note  also,  that  the 
computed  Stanton  number  downstream  of  transition  is  closer  to  measurements 
with  cross  diffusion  included.  Equations  (5.3)  and  (5.4)  indicate  a  transition 
strip  with  fcfl/6,  =  3.0  and  As/6,  =  152.  Numerical  experimentation  shows 
that  transition  can  be  triggered  with  E/,/6,  =  1.0  and  As/6,  =  32. 
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Figure  5.5:  Flow  past  a  prolate  spheroid;  =  20  m/s ec;  with  roughness  strip; 
- New  k-u>  model; - k-u  model  without  cross  diffusion. 
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Figure  5.6:  Flow  past  a  prolate  spheroid;  Uoo  —  30  m/sec;  no  roughness  strip; 
- New  k-ui  model; - k-ui  model  without  cross  diffusion. 
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Figure  5.7:  Flow  past  a  prolate  spheroid;  Uoo  =  30  m/sec;  with  roughness  strip; 
- New  k-u  model; - k-u/  model  without  cross  diffusion. 


x  (  m  ) 


Figure  5.9:  Mach  20  flow  past  a  5°  half-angle  cone; - New  k-w  model; - k-u> 

model  without  cross  diffusion. 
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Chapter  6 

Summary  and  Conclusions 


In  the  context  of  turbulent  boundary  layers,  this  report  shows  that  the  new 
low-Reynolds-number  version  of  the  model  retains  most  of  the  best  features  of 
the  baseline,  high- Reynolds- number  version  of  the  model.  Most  notably,  the 
model  is  just  as  accurate  as  the  baseline  model  for  flows  with  adverse  pressure 
gradient.  This  guarantees  that  the  model  approaches  the  correct  post-transition, 
asymptotic  state. 

The  cross-diffusion  term  is  very  effective  in  eliminating  solution  sensitivity 
to  the  freestream  value  of  ui  for  turbulent  boundary  layers,  free  shear  flows  and 
transitional  boundary  layers.  However,  this  has  been  accomplished  at  the  cost 
of  a  noticeable  increase  in  the  stiffness  of  the  equations  of  motion. 

While  the  new  transition-triggering  method  has  proved  to  be  very  effective, 
and  the  cross-diffusion  term  eliminates  transition  sensitivity  to  the  initial  u> 
profile,  computed  transition  width  is  actually  reduced  somewhat  with  cross 
diffusion  included.  Equation  (5.3)  is  an  upper  bound  for  the  roughness  height 
required  to  trigger  transition  for  all  of  the  cases  considered.  In  many  cases,  a 
smaller  roughness  height  is  sufficient,  and  the  user  should  try  different  roughness 
heights,  if  possible,  to  determine  the  optimum  height  for  a  given  application. 
Equation  (5.4)  is  less  certain.  The  formula  also  expresses  an  upper  bound  that 
covers  all  of  the  cases  considered.  The  actual  values  used  permits  the  roughness 
strip  to  persist  for  at  least  three  streamwise  finite-difference  cells. 
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Appendix  A 


Compressible  Law  of  the 
Wall 


While  the  cross-diffusion  term  is  very  effective  in  eliminating  the  k-ui  model’s 
sensitivity  to  the  freestream  value  of  w,  it  has  a  negative  impact  on  predictions 
for  compressible  boundary  layers.  Figure  A.l  compares  computed  and  measured 
velocity  profiles  for  Mach  4.5  and  Mach  10.3  flat-plate  boundary  layers.  The 
computations  have  been  done  with  ail  closure  coefficients  chosen  to  insure  that 
the  law  of  the  wall  holds  in  the  limit  of  incompressible  flow,  viz, 

u+  =  —£ny+  +  B  (A.l) 

K 

with  k  =  0.41  and  B  —  5.0.  As  shown,  when  cross  diffusion  is  omitted  (trd  =  0), 
the  law  of  the  wall  (in  terms  of  Van  Driest26  variables)  is  reproduced,  and 
close  agreement  with  measured  velocities  is  obtained.  By  contrast,  when  <r<j 
is  0.5  (with  <xd/<r  =  1),  significant  differences  between  computed  and  measured 
velocities  appear  at  both  Mach  numbers.  The  figures  indicate  smaller  differences 
for  ad  —  0.3  (with  <rd/a  —  1/2),  and  this  is  the  value  used  in  all  computations 
of  Chapters  3  through  5. 

To  explain  the  adverse  effect  of  cross  diffusion  on  the  compressible  law  of 
the  wall,  this  appendix  examines  k-w  model  predicted  compressible  log-layer 
structure.  The  results  are  particularly  illuminating  and  clearly  demonstrate 
why  using  too  large  a  ratio  of  <rd  to  <r  in  the  cross-diffusion  term  adversely 
affects  compressible  boundary-layer  predictions.  Aside  from  inclusion  of  the 
cross-diffusion  term,  the  analysis  is  the  same  as  that  presented  by  Wilcox.29 

The  log  layer  is  the  region  sufficiently  close  to  the  solid  boundary  for  neglect 
of  convective  terms  and  far  enough  distant  for  molecular  diffusion  terms  to  be 
dropped,  In  the  log  layer,  the  equations  of  motion  based  on  the  k-ui  model 
simplify  to  the  following. 
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du  , 

^  =  P^T 

(A.2) 

d  \CPT  1  2  .1 

(A.3) 

a*~r  pt  i -  P'puik  =  o 

<tyj  Vdy/ 

(A.4) 

d  f  dwl  p  dk  du)  /du\2  „  ■> 

(A.5) 

p7'  =  p^Tu, 

(A.6) 

The  quantity  uT  is  friction  velocity  defined  as  \frw/pw  where  tw  is  surface 
shear  stress  and  pw  is  density  at  the  surface.  Also,  Tw  is  surface  temperature, 
qw  is  surface  heat  flux  and  Cp  is  specific  heat  at  constant  pressure.  Finally,  y  is 
distance  from  the  surface. 

We  now  change  independent  variables  from  y  to  u.  Consequently,  derivatives 
transform  according  to 

d  du  d  id 

(A.7) 

With  this  change  of  dependent  variables,  we  replace  Equations  (A.3)-(A.5) 
by  the  following. 

d  \CPT  ,  1,  1  qw 

du  \_rrT  2  J  pw  u; 

(A.8) 

••ft 

du2  p2u< 

(A-9) 

d2w  (id  dk  du)  u  (dp^ku)  n 

du 2  it  du  du  ^  °  k  pi,uj 

(A.10) 

Integrating  Equation  (A. 8)  yields  the  temperature,  and  hence  the  density,  as 
a  function  of  velocity  and  Mach  number  based  on  friction  velocity,  MT  =  uT/aw . 


Tw  ^  _1  (7  l^Mr  2(„r)  +  ^«?  («,)+"’  („*)  (A  U) 


Next,  we  assume  a  solution  of  the  form: 

pk  = 


(A. 12) 
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where  T  is  a  constant  to  be  determined.  Substituting  Equations  (A. 11)  and 
(A. 12)  into  Equation  (A.9),  and  noting  that  Mf  =  2VM2,  leads  to  the  following 
cubic  equation  for  T. 


P*  [1  +  (7  -  \)PrT(r* A/T2r]  f2  =  1  (A. 13) 

As  can  easily  be  verified,  when  A, '**  g  1,  the  asymptotic  solution  for  T  is 


— w- — J"'+ 


Finally,  in  terms  of  T,  Equation  (A.  10)  simplifies  to 

_d2W  <Td  dp  ,  r__  ^21  PU 

<T  - — - (-  [o  —  pi  J  - r— 

du 2  p  du  1  pwui  r 


=  0 


(A.14) 


(A. 15) 


Combining  Equations  (A. 11)  and  (A. 12)  yields  the  density  as  a  function  of 
velocity  and  I\ 


Pw 

p 


1  +  (7  -  l)PrT<rTM? 


(A -16) 


Equation  (A.  16)  assumes  a  more  compact  form  if  we  introduce  the  freestream 
velocity,  U0 o  ■  A  bit  more  algebra  yields 


where 


pw _ 1  +  By  -  A'2v 2 

p  1  +  (7  -  \)Ptt<t*YM3 

V  =  U/Uca 

A2  =  ^PttMKToo/T*,) 
U  =  -PrTqw  Uoo/(CpTw  tw  ) 


(A. 17) 


(A.  18) 


Using  Equations  (A.14),  (A. 17)  and  (A.  18),  and  retaining  terms  up  to  0{M 2), 
Equation  (A.  15)  assumes  the  following  form, 


<Puj  ad 

dv3  <7 


where  the  constant  Ku  is  defined  by 


'  B-2A3v  du  1 

1  +  Bv  -  A3v3  dv  J 

[ l  +  Bv-A3v 2. 

K  =  kj  - 


(7-1  )PrT(Za-P/0- 
2  a 


+ 


(A.  19) 


(A. 20) 
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and  k  is  Karman’s  constant.  Because  Uoo/uT  »  1,  we  can  use  the  WKB  method 
to  solve  Equation  (A. 19).  Noting  that  ui  decreases  as  u/U &  increases,  the 
asymptotic  solution  for  u  is 

ui  ~  C  [l  +  Bv  —  A2v2]^1  2<7</a^4exp [—  K^u'/ur]  (A. 21) 


where  C  is  a  constant  of  integration  and  u*  is  defined  by 

«»  _  iol„-i  (  VA2v  -  B  \ 
Uoo  -  A  8m  \\/B2  +  4A2) 


(A. 22) 


Combining  Equations  (A.2),  (A. 12)  and  (A. 21),  we  can  relate  velocity  and 
distance  from  the  surface. 

J  [1  +  Bv  -  A2v2](2",‘/<'~1)/4exp[A'u,u*/uI.]  dv  ~  (A.23) 

We  integrate  by  parts  to  generate  the  asymptotic  expansion  of  the  integral  in 
Equation  (A.23)  as  Uoo/uT  — *  oo.  Hence, 

[1  +  Bv  -  A2v2f+2ai,MI\xp[Kuu-/ur}  ~  (A. 24) 

i  Ut 


Finally,  we  set  the  constant  of  integration  C  =  Tv%/(Kuuw).  Taking  the  natural 
log  of  Equation  (A. 24),  we  conclude  that 


«*  1  , 
- — 

UT  Ku 


+  BU 


(A  .25) 


The  quantity  Bw  is  the  effective  “constant”  in  the  law  of  the  wall  defined  by 


Bu  —  B  + 


PwJ 


(A.26) 


where  B  is  a  true  constant. 

As  discussed  by  Wilcox,29  the  difference  between  Ku  and  the  Karman  con¬ 
stant  is  of  no  great  consequence.  However,  the  variation  of  Bw  with  the  density 
ratio  has  a  large  effect  on  both  skin  friction  and  the  predicted  law  of  the  wall. 
As  clearly  demonstrated  by  Wilcox,13  the  magnitude  of  the  power  to  which 
p/ pw  is  raised  determines  the  degree  of  distortion  of  the  velocity  profile.  For 
example,  the  Standard  k-(  model  is  equivalent  to  a  k-ui  model  with  <Jd  =  2<r. 
This  corresponds  to  an  exponent  of  5/4,  while  the  unmodified  k-w  model  has 
an  exponent  of  1/4.  As  shown  in  Figure  A.l,  there  is  virtually  no  distortion  of 
the  compressible  law  of  the  wall  when  &&  —  0.  By  contrast,  for  the  a  a  —  0.3 
and  0d  =  0.5  cases  we  have  used  <r  =  0.6  and  a  —  0.5,  respectively,  so  that  the 
corresponding  exponents  are  1/2  and  3/4,  respectively.  The  figure  shows  in¬ 
creasing  distortion  of  the  compressible  law  of  the  wall  as  increases.  Although 
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not  shown,  the  k-t  model  departs  even  farther  from  the  compressible  law  of  the 
wall. 

As  a  final  comment,  if  we  had  written  the  cross-diffusion  term  as 

Cross  diffusion  =  (A. 27) 

ui  axj  axj 

the  exponent  in  Equation  (A. 26)  would  be  reduced  to  1/4,  independent  of  the 
value  of  (Td-  This  is  true  only  for  the  constant-pressure  case  since  pk  is  ap¬ 
proximately  constant  in  the  log  layer.  The  same  statement  cannot  be  made 
for  boundary  layers  with  nonzero  pressure  gradient,  however.  A  shortage  of 
Contract  funds  precluded  testing  this  alternative. 
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